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ABSTRACT 



The interaction of a sound field and a seamount has been 
studied by physically modeling the ocean surface over 2 
dimensional and 3 dimensional models of Dickens Seamount. 

By using a smooth ocean surface, and one with a surface of 
scaled Rayleigh roughness to model a 35 knot wind, the 
relative contributions are determined for off-axis scattering 
elements, multiple reflection from the ocean surface, and 
diffraction over the crest of the seamount. 

Boundary wave generation over a randomly rough plane sur- 
face is studied experimentally. The ratio of boundary wave 

amplitude to volume wave amplitude is found to be proportional 
3/2 1/2 

to (frequency) and (range) and the ratio of the empiri- 
cal scattering parameter to the rms height of the roughness 
elements is found to be approximately 0.3. The spatial corre- 
lation length of the randomly rough surface takes the place 
of the center-to-center separation of the hemispherical bosses 
used in Tolstoy's theoretical treatment. 
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INTRODUCTION 



I . 

In previous work the interaction effect of an acoustic 
wave with an underwater barrier such as a seamount has been 
considered in terms of forward scatter at the upslope and 
diffraction over the crest. This work was grounded on the 
theoretical developments of Biot and Tolstoy [Ref. 1] which 
provided a solution to the problem of sound energy diffracted 
by a rigid, infinite wedge. Bremhorst [Ref. 2] was able to 
demonstrate close agreement between the Biot-Tolstoy theory 
and experimental measurements for a 90 degree wedge used as 
a diffracting barrier. 

Spaulding [Ref. 3] extended this approach to a physical 
scale model of Dickens Seamount located in the Gulf of Alaska. 
In order to validate the laboratory model, Spaulding used the 
ocean data of Ebbeson et al . , [Ref. 4] of the Canadian Defense 
Research Establishment, Pacific (DREP) . Figure 1 supplied by 
DREP shows their 1978 ray patterns for three cases of inter- 
action of the sound field with Dickens Seamount. Spaulding 
addressed only case A where the rays appear to be completely 
blocked. Figure 2 shows that for case A there was a dis- 
crepancy of about 15 dB between theory and what was observed 
at sea. Concurrently, Medwin and Spaulding [Ref. 5] were 
able to extend the Biot-Tolstoy theory, show how it could be 
used for a finite wedge, and define the Diffraction Strength 
of a wedge-like crest as shown in Figure 3. However, in the 
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FIGURE i. Ray traces for track b. Source ranges are 

79 Km (A) , 99 Km (B) , and 119 Km (C) . Source 
depth = 18 m, receiver depths = 329 m, 633 m. 
The source angles are ±15° in 1° increments 
(from Ref. 4) 
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FIGURE 2. Measured propagation loss for track 6 with the source at 18 m and 
receiver at 329 m. The smooth line is the FACT prediction based 
on a flat bottom. The arrow indicates the seamount range (from Ref. 4) 
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Diffraction strength, Definition (from Ref. 5) 
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work, Spaulding did not consider the effects of the inter- 
action of the sound energy with the ocean surface above it. 

In May 1978 Tolstoy [Ref. 6] presented a new theory of 
acoustic scatter from a slightly rough surface at near graz- 
ing incidence which incorporated the boundary condition for 
a rough surface developed by Biot [Ref. 7] . This theory 
represented a significant departure from more traditional 
approaches to the problem of acoustic scatter for three pri- 
mary reasons: First, the theory is first order in the rough- 

ness parameter, as opposed to more conventional approaches 
which consider the roughness parameter a second order effect: 
Second, the theory includes multiple scatter and diffraction 
which is generally ignored by other theories: Third, Tolstoy's 

theory does not make the Kirchhoff Assumption which is par- 
ticularly untenable at low grazing angles. 

The surprising result of Tolstoy's theory is the predic- 
tion of the existence of what he called a "boundary wave" in 
the fluid above a slightly rough surface generated at near- 
grazing incidence by coherent multiple Rayleigh scatter. 

1/2 

Since the r power law Tolstoy derived for the range de- 
pendence of the ratio of the boundary wave to the volume 
wave (the volume wave is that wave which exists in the fluid 
well above the scattering surface) was that of cylindrical 
spreading instead of the spherical spreading which the volume 
wave suffers; the theory also predicted that as the range 
increased, the boundary wave would become stronger than the 
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volume wave at some range. This effect was first observed 
by Bailie [Ref. 8] for the case of scattering by closely 
packed hemispherical bosses on a rigid planar surface. The 
agreement between Bailie's results and Tolstoy's theory was 
quite remarkable. 

The next author to consider the problem was Hollis [Ref. 

9] . He addressed himself to a combination of the wedge prob- 
lem and the coherent scatter problem by applying the hemis- 
pherical bosses used by Bailie to a 14 degree rigid wedge. 

The slope of 14 degrees was chosen because that corresponds 
to the average slope of Dickens Seamount. Hollis' results 
closely matched those of Bailie and again verified the pre- 
dictions of Tolstoy by demonstrating the existence of the 
boundary wave. His results also showed for the highly 
idealized case of the 14 degree wedge roughened by hemispheri- 
cal bosses, that the energy contained in the boundary wave as 
it diffracted over the crest of the wedge could be larger than 
the diffraction of the volume wave. 

The purpose of this work is to consider the effects of 
including an appropriately roughened ocean surface above the 
model of Dickens Seamount for the geometry used by Spaulding 
and to investigate the possibility of extending Tolstoy's 
boundary wave theory to a randomly rough seamount surface. 
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II. THEORY 



A. DIFFRACTION OVER A WEDGE 

In general, diffraction occurs whenever an acoustic wave- 
front encounters an obstacle of any sort. Of course, whether 
the diffraction effects are important or if the phenomenon 
being observed can be adequately represented by the approxi- 
mations provided in ray theory, is a function of the specifics 
of the problem under study. For the case of a source and 
receiver separated by an infinite, rigid wedge with no direct 
path between the two, diffraction provides the only mechanism 
by which acoustic energy emanating from the source can reach 
the receiver. 

Consider the case of an infinite wedge represented in 
Figure 4. The wedge is made up of two perfectly rigid plates, 
which intersect at a Crestline. This crest will provide a 
convenient geometrical orientation and is thus designated as 
the z axis of a cylindrical coordinate system. The region 
in space not occupied by the wedge is filled with a homoge- 
neous, compressible fluid of density p where sound travels 
with velocity C. In such a system the displacement potential, 
<5, is described by 






1 

r 3r 



2 







( 1 ) 



Harmonic solutions of this equation have the form 
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FIGURE 4 . Geometry of Biot-Tolstoy theory for diffraction 
by a wedge. Wedge angle 9 is measured in the 
fluid (from Ref. 5) ^ 
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and 



(3) 



(1 2 ) 

H ' are Bessels functions of the first 
kind 



The rigid plates that make up the wedge require 



=0 at 9 = 0 and 9 = 9 

9 0 w 



(4) 



Applying this boundary condition to Equation (2) implies 



= cos V 9H 
n V 



(5) 



where 



mr 



n 



— — , n— 0/1,2/... 



w 



( 6 ) 



Since at r = 0 the imaginary part of the Hankel function is 
infinite, it is further required that only the real part of 
will be retained. 

V 

Now assume a point source S (volume/time) which starts 
to flow uniformly and instantaneously at t = 0 , so that at 
range R after time t 



l± = (-S/4ttR) 1 (t -R/C) (7) 

3t 
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Equation (5) is symmetric about the z axis. The nomal 
coordinate method of Biot-Tolstoy is followed (see Ref. 1 
for a detailed derivation of the method) . Thus, the solution 
to Equation (1) becomes 

3 ^ $ C r 

— o = - ) COS V 9 COS V e^I (8) 

^ n 0 n 

dt w n 

for ct ^ z, where 

00 

(kr) (kr^) [k (C^t^-z ] k dk (9) 

Using the transform detailed in Ref. 1, Equation (8) reduces 
to a form that is more easily visualized. Let 



If/ N 2 ^ 2,1/2 

c‘(r -ro> +21" 



If/ ^ \ 2 ^ 2,1/2 

^[(r+r^) + z ] / 



( 10 ) 



t^ is the time of travel of a pulse of sound from source to 
receiver by direct path. Tq is the time of travel from the 
source to the wedge crest then to the receiver. Thus, Xq is 
the time of travel of diffracted energy. This separates time, 
and therefore Equation (8), into three distinct regimes; 
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^0 
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Since it is diffracted energy being considered, this implies 
solutions must be valid for Tq < t. For this case 



X -1 . n 

I = : — r — sin V e 

n ^ sinh y n 



( 11 ) 



where 



y = 



2 2,2 2 2 , 

, c t - (r + r » + z ) 

cosh"-^[ ^ 



2 rr 



( 12 ) 



0 



Substituting the value for back into Equation (8) 



3 Q 



-v„y 



3t 



2 TT0 r r,^ sinh y 



y cos V 9_ cos V 9 sin v ire (13) 

- n U ^ ^ 



n 



n 



where < t and y is given by Equation (12) . Expressing the 
trigonometric functions as exponentials 




^ exp(-7ry/0^) 

4 tt9 r r_ sinh y 
w 0 



sin (tt/9, ) (tv±9±9„) 
w u 

1-2 exp (-Tvy/9^) cos (^/9^) 
• (Tr±9±9Q)+exp(2iTy/9^) 



(14) 



Since the acoustic pressure, p, and the displacement poten- 
tial, ■$, are related by 



M T 

3t 

The acoustic pressure due to the diffracted wave in the 
shadow zone of the wedge is 
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_pC exp(-7:y/9^) 

4 7:0,, r r- sinh y 

sin(:r/0 ) (7rt9±9 ) 

w U 

l-2exp(-:ry/9 )cos(tt/9, ) {t:±0 + 9 -) +exp (-2Try/0 ) 
w w 0 w 



( 16 ) 



where 



y = cosh 



-1 



^ 2^2 . 2 ^ 2 . 2 , 

C t - (r + Tq + z ) 

2 r r„ 



(17) 



and 



9^ = fluid region above the wedge 

0Q = angle between wedge and source 

0 = angle between wedge and receiver measured 

from source side of the wedge in the fluid 
above the wedge 

r = distance from source to crest of the wedge 

z = distance along the crest of the wedge where 
z = 0 represents the least time path from 
source to crest to receiver. 



As pointed out by Spaulding [Ref. 10] the diffraction 
phenomenon can be described qualitatively as follows: The 

acoustic wave produced by a point source at a distance r from 
the crest of the wedge expands spherically. The wavefront 
first encounters the wedge crest at the point of least time 
travel path at time t = r/C and z = 0. The wave will then 
continue to interact with the crest of the wedge at increasingly 
later time and increasing Izj. In effect, the crest of the 
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wedge acts as a continuous, time shaded, line source. Medwin 
[Ref. 11] has computed the frequency transform of Equation 
(16) . It was this transform that Spaulding used to calcu- 
late theoretical values to be compared with experimental 
data for three models of Dickens Seamount of increasing 
complexi ty . 

B . BOUNDARY WAVE 

I. Tolstoy has presented an application of a theory 
developed by M. A. Biot [Ref. 12] foinnulating sound scatter from 
slightly rough surfaces at near-grazing incidence. Tolstoy's 
theory [Ref. 6] yields solutions for the coherent scatter of 
transient spherical waves by rough planes where the size of 
the roughness is small compared to the wave length of the 
signal. Biot had previously demonstrated that if the spacing 
between hemispherical bosses was small compared to a wave 
length, they could be replaced with continuous distributions 
of monopole and dipole radiators [Ref. 13]. 

To extend this formulation, Tolstoy assumed: 



d 

r 



< 



h 

X 



< < 



1 



(18) 



where 



d = diameter of hemispherical boss 
h = center to center spacing of bosses 



A 



acoustic wavelength 
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By applying the boundary condition formulated by Biot 
[Ref. 14], Tolstoy was able to show that for near-grazing 
incidence a boundary wave would be developed and that this 
boundary wave would decrease exponentially as the receiver 
was moved away from the plane of the hemispherical bosses. 
The temporal solution describing the pressure of the boun- 
dary was then transformed to the frequency domain using 
Fourier Transform techniques giving the pressure in the 
boundary wave as a function of frequency. The details of 
this development are contained in Ref. 6. 

Thus, Tolstoy reported in 1979 that the amplitude of the 
postulated boundary wave as a function of frequency is 

^BWA ^ -^[jQ(kr) + (kr) ] ^^^u^exp [-e (z+Zq) k^] (19) 

where 

£ = scattering parameter and is proportional to 

the volume of scatters per unit area 

k = wave number 

z = height of receiver above the plane 
Zq = height of source above the plane. 

For the far field case where kr >> 1, the entire scattered 
field can be written as 

P (f) = k^{e^ exp[-2ek^ (z+z ) ] 

s 2tt kr 0 

W-, 9 1/2 

+ 2 (^) ■^/'^6 sin (kr-^) exp [-ek (z+z^]} (20) 
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where, for closely packed hemispherical bosses with center 
to center distance, h, of 2 ram in air. 



t, — 



8 . 88 X 10"^ m 



k = 



2Trf ^-1 
345 ^ 



r = range from source to receiver in meters 



f = frequency in Hertz 



2 - = height of source in meters 



height of receiver in meters 
2 + z. 



9 = arc tan(- 



'0 



For the case of the source and receiver in the plane of the 
roughness then, z = z^ = 0 and Equations (19) and (20) reduce 
to 



?3(f) 






^ ,,2 ^ 1/2 

2tt ^ ^kr-' 



.3/2 



/2 



TTr 



( 21 ) 



If the source is an impulse, then the acoustic pressure 
generated by that source is 



'‘5 • 



( 22 ) 



Transforming to the frequency domain 
^VWA ^ ^ ^ Jnr 



(23) 
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If Equation (21) is divided by Equation (23), an expression 
of the ratio of the boundary wave to the volume wave (spher- 
ically spreading wave) is found which represents the frac- 
tional change in wave amplitude in the plane of the hemispherical 
bosses caused by those bosses 



^BWA 

^VWA 



e (2Trr) 



1/2 






(24) 



It is this expression that was compared against experimental 
data by Bailie [Ref. 15], Medwin, et al., [Ref. 16], and 
Hollis [Ref. 17] and found to be in excellent agreement for 
kh < 1. 
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Ill . RESEARCH FACILITIES AND EQUIPMENT 



A. OCEAN PHYSICS LABORATORY AND ANECHOIC CHAMBER 

The portion of the investigation dealing with the scale 
model of the Dickens Seamount was conducted in the Ocean 
Physics Laboratory located at the Naval Postgraduate School. 
The excellent signal to noise ratio provided by the source/ 
receiver selection made it unnecessary to use extraordinary 
quieting procedures, other than a relatively large number of 
signal averages prior to signal processing. The signal to 
noise ratio for this part of the experiment was never less 
than 20 dB from 20 kHz to 100 kHz for all of the work done in 
the laboratory itself. This fact allowed the investigation 
to proceed without requiring the seamount model to be moved. 
As the model is quite bulky and massive, this alone greatly 
facilitated the procedure. 

On the other hand, that portion of the investigation 
dealing with forward scatter from a randomly rough plane sur- 
face was conducted in the Anechoic Chamber located next to 
the Ocean Physics Laboratory. This was a result of several 
factors all dictating a change in the source selection and 
the resultant decrease in available signal strength. As 
reported by Hollis [Ref. 18] receiver 4132 (with no source) 
in the Anechoic Chamber compared to the same receiver blocked 
by two sheets of 1/4 inch aluminum showed noise level differ- 
ences of 0.1 dB or less over a frequency range of 2-30 kHz, 
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suggesting that the noise was electronic rather than 
acoustic . 

B. DATA ACQUISITION AND PROCESSING EQUIPMENT 

Data acquisition and processing was accomplished for both 
portions of the investigation by a digital computer system 
composed of four major components which were interfaced to 
provide high speed analog to digital conversion, data pro- 
cessing, data storage and output. The design, itself, was 
developed by the Special Projects Section, Naval Air Develop- 
ment Center in cooperation with Pinkerton Computer Consultants, 
Inc. of Warminster, Pennsylvania. The primary components of 
this system are as follows. 

1 . Interdata Model 70 Computer 

This minicomputer has a 64 thousand byte core memory 
and is programmable in FORTRAN or BASIC. In addition to core 
memory, information can be stored on floppy disk or cassette 
tape for later processing. 

2 . Phoenix Analog to Digital Converter, Model ADC 712 
The Model ADC 712, is a high speed, high accuracy 

analog voltage to digital converter which can encode a 20 
volt, peak to peak, input signal into 12 binary bits, with a 
resultant resolution of one part in 4,095 at a maximum rate of 
2 microseconds per conversion. This fast settling and digiti- 
zation by encoding process allows a typical computing through- 
put rate of 476,190 channels per second, including settling 
time. The sampling frequency is provided by a General Radio 
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Decade Oscillator, via a locally developed sampling circuit. 
The stability of the oscillator is 0.001% of indicated value. 



3 . Texas Instruments Silent Electronic Data Terminal, 

Mode 1 733 

The TI 733 consists of a keyboard used as a printer, 
a programming input/output control device, and a transmit/ 
receive device interfaced with a peripheral, dual disk, drive 
unit. The system provides rapid, highly accurate processing 
and was primarily used for both time and frequency domain 
analysis. The frequency domain analysis was effected by 
standard Fast Fourier Transform techniques . 

4 . Qrbis Model 76 Diskette Drive 

The Orbis Model 76 is a compact, portable, dual drive, 
direct access, 256 thousand byte floppy disc data storage 
device. The floppy magnetic discs provided the capability 
of long term storage of vast quantities of data for later 
analysis . 

C. EQUIPMENT LIST 

The following is a list of the equipment used during the 
course of the investigation. Also, the abbreviated name of 
each piece of equipment is listed as it will be referred to 
throughout the text. 



28 



TABLE I 





EQUIPMENT LIST 


Description 


Abbreviation 


Scope 


Tecktronic Type 551 Dual-Beam 
Oscilloscope 



Frequency Synthesizer General Radio 1312 Decade Oscillator 



Timing Simulator 


Interface Technology Timing Simu- 
lator/Word Generator, Model RS-648 


Frequency Counter 


Hewlett Packard Model 522 3L 
Electronic Counter 


Filter HP 


Krohn-Hite Frequency Filter, 
Model 3342 


Filter BP 


Krohn-Hite Frequency Filter, 
Model 3322 


PAR Amplifier 


Princeton Applied Research Ampli- 
fier, Model 113 


MIC Power Supply 


Bruel and Kjaer Microphone Power 
Supply, Model 2804 


ARB 


Wavetek Arbitrary Waveform 
Generator, Model 175 


MIC Preamp 


Bruel and Kjaer Model 2619 


Power Supplies 


Hewlett Packard Model 721A 


Power Amplifier 


Kilowatt Amplifier Model L2 
Instruments, Inc, 


Oscilloscope 


Hewlett Packard Model 14 OA 
Oscilloscope 


Nicolett 


NSC 660B Dual Channel FFT 
Analyzer 
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FIGURE 5. Ocean physics laboratory equipment 
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FIGURE 6. Anechoic chamber equipment 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 



A. SEAMOUNT 

The physical models used for this investigation were the 
same ones Spaulding used in 1979 [Ref. 19] . They were of 
differing levels of complexity. The simpler of the two, the 
2D (PE Contour) , was constructed using a technique similar to 
that used in the fabrication of aircraft wings [Ref. 20] . 

It conformed to the profile of Dickens Seamount along track 
6 of the ocean traversed by Ebbeson [Ref. 21] . 

This model had the profile of the seamount along track 6 
but had none of the variation of the real seamount in a direc- 
tion perpendicular to track 6. As such, this model repre- 
sented a physical manifestation of the 2D, parabolic equation 
computer model used by Jensen [Ref. 22] . 

The more complex of the two models, the 3 dimensional 
scale model of Dickens Seamount, was also constructed for 
Spaulding in 1979 [Ref. 23]. This model was constructed by 
Modelmakers, Inc., of San Francisco to conform to the bathy- 
metric data supplied by Ebbeson. The seamount model was con- 
structed by laminating 1/8 in layers of particle board cut to 
represent the shape of the real Dickens Seamount as accurately 
as possible at the appropriate depth. The result, to this stage 
in construction, was what would appear to be a contour 
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map of the seamount if it were viewed from overhead. Once 
all layers were completed and assembled in the appropriate 
sequence, the entire structure was covered by a thin layer 
of plaster to smooth the transition from layer to layer and 
to add surface density to the model. The final result was 
a three dimensional, acoustically rigid, scale model that 
accurately reflected the known bathymetry of Dickens Seamount 
to a resolution of what would correspond to 25 m at sea. The 
scale used was 5 inches to represent 1 km or 1:7874. 

B. RANDOMLY ROUGH PLANE SURFACE 

Three randomly rough surfaces were constructed of aluminum 
plates covered with pea gravel. The plates were 48 inches 
by 60.75 inches, 48 inches by 48 inches, and 19.7 inches by 
19.7 inches. The 48 inch by 60.75 inch surface was used to 
model a randomly rough ocean above the model of Dickens 
Seamount. The 48 inch by 48 inch surface was used to extend 
the work done by Bailie [Ref. 22] for a surface of hemispheri- 
cal bosses to the case of a randomly rough plane surface. 

The 19.7 inch by 19.7 inch surface fitted with a height 
measuring micrometer and track was used to investigate the 
statistics of these surfaces. 

The construction procedures were very straightforward. 

The surface involved was covered with melted paraffin to a 
depth of approximately 0.5 cm. Paraffin was chosen because 
of its ease in application and low melting point. As Hollis 
pointed out [Ref. 24], the paraffin serves much better than 
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other materials because of the ease with which one can 
eliminate bubbles in the substance and thus avoid the problem 
of bubble resonances. After the paraffin had cooled a short 
time, pea gravel was added until a relatively thick layer that 
completely covered the paraffin was achieved. Then, using a 
hand-held hot air blower, the surface was reheated until the 
paraffin began to melt. As the paraffin melted, the gravel 
would sink down into it. Thus after the paraffin had cooled 
a second time, the gravel was securely imbedded in it and 
therefore attached to the aluminum plate. After the entire 
plate was cooled a second time, touch up additions were made 
to any area in the gravel surface where a significant amount 
of paraffin was visible. When the procedure was complete, 
there was no area on any of the three plates where the pre- 
sence of the gravel was not the dominant feature of the 
surface . 

As the purpose of each of the plates was different, so 
were the details of each plate before the gravel was applied. 
The plate to be used to model the ocean surface was 48 in by 
60.75 in. The gravel surface itself covered an area of 45 in 
by 55 in on this plate. After the gravel was applied, the 
model ocean was then bolted to a supporting frame to minimize 
curvature when suspended and reduce the effort required to 
move it. Then the plate was suspended at a height appro- 
priate to the scaling of the seamount using a chain fall. A 
detailed look at the specific appearance of this model ocean 
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while in place above the seamount is contained in Figures 
7 and 8 . 

The plate to be used for boundary wave investigation was 
48 inches by 48 inches. In the center of the plate a hole 1.5 
inches in diameter was drilled to admit the source. Then 
holes 0.32 cm in diameter were drilled at ranges of 10 cm, 

20 cm, 30 cm, and 40 cm from the source hole along one radius. 
This drilling procedure was repeated along radii that were 
separated by 15 degrees in a complete circle around the source 
hole. Thus there were a total of 96 holes of 0.32 cm diameter 
in the plate surrounding the source hole. The holes were 
plugged from the under side of the plate and gravel applied 
to the upper side. The details of the appearance of this 
plate with source and receiver in place are contained in 
Figures 9 and 10 . 

The sample plate used for statistical study measured 50 cm 
by 50 cm. After the gravel was applied to this plate, it was 
installed in a locally designed apparatus to measure the 
surface characteristics. This device was constructed to allow 
the placement of a depth micrometer anywhere in the x-y plane. 
A 5 cm by 5 cm area was then chosen as representative of the 
entire plate and the depth measured at intervals of 1 mm. 

These data were formed into a square matrix containing 2600 
elements from which the probability density of the heights and 
slopes were calculated. As can be seen in Figures 11 and 12, 
the probability density functions for both the heights and 
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FIGURE 7. Seamount with model ocean in place 







Upslope side of seamount model with model 
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tIGURE 9. Randomly rough plane surface, gravel side 
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FIGURE 10. Randomly rough plane surface, control side 
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SURFACE SLOFE FROBABILITY DENSITY FUNCTION 



slopes are very nearly Gaussian. The rrns height was calcu- 
lated to be 2.23 mm and the rms slope was calculated to be 
1.15. The correlation function was calculated by row and by 
column. This was done for each row/column; then the average 
taken. Figure 13 shows the plot of the results of the aver- 
age for both rows and columns. If a correlation length is 
defined as 

C(L) = e~^ (25) 

then the spatial correlation length of this randomly rough 
surface, and therefore all three plates, is 



L = 4 . 8 mm 



(26) 



As can be seen from Figure 13, over this distance there is 
virtually no difference between correlation lengths as calcu- 
lated by row or by column. 

To validate the use of this surface to model the ocean 
surface we assume that the 35 knot wind during the ocean 
experiment generated a Pierson-Moskovitz wind/wave spectrum 
[Ref. 25]. That semi-empirical spectrum is 



$ (n) 




exp[-3(-j^) ] 



where 

a = 8. 1 X 10 

3 = 0.74 



(27) 



42 



o □ 



(X3 




a 






1 

1 


i t 1 




o 


00 


<S) 




C\J 


o 




d 


O 


d 


O 


d 



C\J 

O 

1 



43 



p By Column 

AVERAGE SPATIAL CORRELATION FUNCTION 



^19 5 ~ Wind speed at 19.5 mm above surface 

g = acceleration due to gravity. 

Then the rms height of the waves is given by 

» 1/2 

a = [ / 4(fi)dn] (28) 

0 

ct 

= ^ = 1.74 m for 35 knot wind 

4sg 

Since the geometry in the laboratory with the seamount 
model was the same as the geometry at sea during the CW 
experiments, modeling of the Rayleigh roughness parameter 
requires only that ka be the same for the two situations. 



[ka 



ocean 



[ka] 



lab 



lab 



[fa] 



ocean 

lab 



for f = 50 Hz 

ocean 



lab 



= 18.1 kHz 



(29) 



for f 

ocean 



500 Hz f^^j^ = 181.3 kHz 



It is assumed that the phase shift at the pressure release 
ocean surface does not have to be modeled in the laboratory 
surface because of the randomness of the surface scattered 
signal . 
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C. SOURCE/RECEIVER SELECTION 
1 . Seamount 

Since the objective of the analysis was to extend the 
work Spaulding has done [Ref. 26] the choice of source/ 
receiver was made to correspond to that made by him. This 
choice was also dictated by the ocean experiment of Ebbeson 
[Ref. 27]. As shown in Figure 14, the bulk of energy [Ref. 28] 
intersected the source side of Dickens Seamount about midslope. 
This energy was contained in two major bundles which apparently 
divided slightly before intersecting the seamount itself. In 
an effort to duplicate this effect as closely as possible an 
extended source was used. That is, a low "Q" 8 cm by 9 cm 
rectangular solid dielectric transducer was selected as the 
source and a 1/2 inch B and K microphone as the receiver. 

These were of the same type used by Spaulding during that 
portion of his procedure [Ref. 29] . The choice of receiver 
was driven by the twin requirements for both small size, and 
thus increased high frequency response, and sufficient sensi- 
tivity. The choice made was thus the 1/2 inch B and K micro- 
phone, its associated preamplifier and power supply. In some 
ways, this choice of source was unfortunate, as it was later 
discovered, since this type of source proved to be unstable 
both in time and in frequency. The variability shown in 
Figures 15 through 18 demonstrate the magnitude of the problem. 
The attempts at finding a solution to this problem will be 
discussed at length in the section dealing with signal 
processing . 
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Randomly Rough Plane Surface 



2 . 

Equation (24) assumed a point source and point re- 
ceiver both on the surface. The selection of a source/receiver 
combination became a tradeoff between size (since ka << 1 
where a is the radius of the source to approximate a point 
source) and signal strength. In view of this tradeoff, the 
B and K Model 4145 one inch diameter microphone was chosen for 
the source and the 1/2 inch B and K Model 4133 with a 3 . 2 mm 
probe tube was chosen for the receiver. 

To increase the radiation parallel to the plate the 
source was mounted in a circular cap 2.5 cm in diameter with 
a hole 1.3 cm in diameter in the side. The cap was filled with 
tacky wax to a level and at an angle such that the existence 
of a standing wave was not possible. The cap itself was part 
of a larger device that permitted vertical adjustments using 
a micrometer and rotational adjustments using a compass rose. 
Figures 9 and 10 show the source, both from the upper and 
lower sides of the randomly rough plate. 

The half inch B and K Model 4133 receiver was attached 
to a 3.2 mm probe tube that was fitted through a micrometer 
dial. It also had a nylon spacer between the barrel of the 
probe tube and the cap that actually fitted to the 1/2 inch 
B and K microphone itself. This spacer provided electrical 
isolation between the aluminum plate and the microphone and 
therefore betv/een the source and receiver. The probe tube 
v;as 7.8 cm long from the end of the tube to the nylon spacer. 



47 



i 

1 

i 



1 




) 




The small diameter of the probe tube ensured that the receiver 
was small compared to a half wave length for the highest fre- 
quency of interest to avoid pressure averaging over the face 
of the microphone. This was a potential problem because of 
the flush mounting of the receiver. 

The power supply for the receiver was battery powered. 
This significantly reduced the electrical noise in the system. 
Since the heating element was not used, thermal noise was 
also significantly reduced. Since the height of both source 
and receiver was set by micrometer, firm control of the rela- 
tive source/receiver height was obtained. This was important 
because, as reported by Bailie [Ref. 30] and Hollis [Ref. 31], 
the boundary wave was sensitive to source height relative to 
the roughness elements. 

D. SIGNAL PROCESSING 
1 . Source Signal 
a. Seamount 

Because of the response of the rectangular solid 
dielectric transducer, a significant amount of experimentation 
was required to find the most advantageous combination of 
polarization voltage, driving voltage, and wave form. For- 
tunately, the Wavetek 175 Arbitrary VJaveform Generator (ARB) 
was designed to provide the flexibility necessary to adequately 
address these complex problems. 

First it was noticed the wave form produced by 
the source was fairly sensitive to the polarization voltage 
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applied. This was especially true of the time necessary to 
insure that the ringing had decayed to a level that could be 
considered negligible. After several different voltages were 
examined, the obvious best choice was +150VDC. This produced 
a signal that was surprisingly free from ringing when con- 
sidered in the context of the low "Q" of the transducer. The 
signal, as viewed without obstruction of the seamount is con- 
tained in the top graph of Figure 19 . 

Second, the amplitude of the driving signal was 
addressed. The problem was considered in the context of 
getting maximum output from the transducer without overdriving 
it. Again, after a significant amount of experimentation, a 
driving voltage of 110 V AC proved to be of shortest duration 
for the polarization voltage applied. 

The third, and probably most difficult, considera- 
tion was the actual wave form to be applied to the transducer. 
At this point in the effort, the tremendous advantages the 
ARB afforded came to the fore. In addition to several pre- 
programmed wave forms stored in Read Only Memory (ROM) , the 
ARB also provides several Random Access Memory (RAM) devices. 
Since the resultant transducer response to the preprogrammed 
waveforms proved unsatisfactory, several waveforms were pro- 
grammed in the RAM. Each RAM consists of 256 discrete blocks 
of memory. The voltage to be applied via the ARB is then 
programmed by the user in each of these 256 blocks. The 
voltage is then entered as a quantized level from 0 to ±127. 
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The actual voltage that is assigned to each quantum level is 
controlled by the user via the amplitude setting provided. 

Thus for an amplitude setting of 1 volt, +127 would corres- 
pond to +1 volt and -127 would correspond to -1 volt. Using 
this approach it was observed that the best output from the 
dielectric transducer was obtained by programming a RAxM to 
apply +127 from the first block, -127 from the second block, 

0 from the remaining 254 blocks, setting the amplitude to 
0.5 volts, applying a frequency of 1000 Hz to the process; 
and, finally, amplifying the results to a value of 110 VAC. 
Therefore the voltage applied to the source was a positive 
impulse of 110 VAC followed by a negative impulse of -110 
VAC. The voltage was actually applied to the source for 7.8 
microseconds. This process was repeated every 25 milliseconds. 
A detailed list of equipment settings is contained in Table II, 
Appendix A. 

b. Randomly Rough Plane Surface 

In order to be used as an acoustic source, the 
one inch B and K microphone must first be biased with a DC 
voltage. Once this has been accomplished, the AC signal is 
imposed on top of the bias voltage. The B and K Model 4145 
one inch microphone is rated by the manufacturer at a voltage 
of 250 V (peak plus bias) [Ref. 32] . Thus with a bias voltage 
of +150 VDC an applied AC signal of 100 VAC represents a 
driving voltage that is as large as possible yet still within 
the manufacturer's specifications. Therefore, this combination 
of DC and AC voltages was used. 
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In an effort to get as much energy in the har- 
monics of the driving frequency as possible, a half triangu- 
lar pulse was selected from those preprogrammed in the 
Wavetek Arbitrary Waveform Generator (ARB) . This pulse was 
driven at a fundamental frequency of 2.5 kHz. This gave a 
frequency spectrum with relatively significant amounts of 
energy at frequencies that greatly exceeded kL 1, where 
L is the spatial correlation length of the gravel surface. 

2 . Received Signal Processing 
a. Seamount 

With the two models of Dickens Seamount in the 
Ocean Physics Laboratory, an increase in the signal to noise 
ratio was achieved by taking 1500 samples of each signal 
being considered. Since, for N samples of a signal in the 
presence of noise 

A(-) = 10 log /N (30) 

n ^ 



Thus, for 1500 averages a gain of 16 dB was achieved in the 
signal to noise ratio. 

A second consideration that needed to be resolved 
before actual data taking began was that of the sampling fre- 
quency. As stated by the uniform sampling theorem, if any 
signal which is band limited is sampled in time, with a uni- 
form intersample period T^ such that 



T 



s 




(31) 
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where 



W = Bandwidth 

the signal is completely specified [Ref. 33]. This leads to 
the result 




s 



But, given reality, if a figure of 3W is used rather than 2W, 
the data collected is much more reliable. Since the highest 
sampling frequency with which the A/D equipment could deal was 
320 kHz, this would allow reliable frequency data up to 
approximately 100 kHz. 

Therefore, using a sampling frequency of 320 kHz, 
1500 averages, and a pulse repetition rate of 25 milliseconds, 
the data from the 2D (PE Contour) and the scale model Seamounts 
were collected and stored on magnetic disk for later processing. 
These recorded signals were of both seamount models without 
ocean, both seamount models with mirror ocean, and 21 runs 
for both models with rough ocean. Since the special correla- 
tion length had been established, L = 0.48 cm, by shifting 
the position of the ocean above the seamount model a statis- 
tically "new" ocean had, in effect, been placed over the model. 
By shifting the ocean that had originally been centered over 
track 6 ± 10 cm in 1 cm increments, 21 statistically indepen- 
dent manifestations of the randomly rough ocean had been 
placed over the model. This allowed statistical methods to 
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be used to arrive at an average pressure as a function of 
both time and frequency. 

It was during the taking of the data with the 
rough ocean in place above the scale model seamount that possi 
ble problems with the source became apparent. While observ- 
ing the plot of pressure versus time a change in level as a 
function of run number was noticed. Since the first peak 
in amplitude of the received signal was entirely due to 
diffraction, the presence of the ocean should have made no 
difference in the amplitude of this first peak. The ampli- 
tude of the first peak as a function of run number is shown 
in Figure 15. The entire procedure for all 21 runs took 
place over a period of approximately four hours. As this 
figure demonstrates, apparently there was a very low frequency 
oscillation in the source strength presumably due to instabili 
ties in the source itself. This immediately raised the ques- 
tion of the frequency stability of the source. 

With the rough ocean in place above the seamount, 
the duration of the first peak is represented by somewhat 
more than 16 data points. Since the algorithm used in the 
Fast Fourier Transform requires the use of the number of data 
points such that the number is an exact multiple of 2, the 
first 16 points of the pressure versus time signal that was 
above background noise was transformed into the frequency 
domain. This led to a frequency resolution of 20 kHz. After 
choosing the case of the 2D (PE Contour) model with the 
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smooth ocean an arbitrary reference to which to normalize the 
frequency data representing the first 16 points of each sig- 
nal that involved an ocean (rough or smooth) was plotted on 
the same graph as the reference. The objective was to identify 
a correction factor to be applied to a particular run based 
on the frequency behavior of the first amount of energy of 
that which was known to involve diffraction only. Three 
examples of this procedure are contained in Figures 16, 17, 
and 18 . 

As can be seen from these three figures, the fre- 
quency response of the source in the neighborhood of 30 kHz 
is relatively stable. However, above 40 kHz the variance from 
run to run was extreme. It was also apparent from the be- 
havior of Spaulding's frequency plots that the far-field of 
the system was not reached until approximately 50 kHz. Based 
on the foregoing, it was concluded that the time data was 
reliable for approximately 30 kHz but not necessarily for the 
higher frequencies. It was concluded that any information 
about the effect of an ocean above a model of the seamount: 
and its effects on the acoustic shadowing by that seamount 
would have to be derived from the pressure versus time data. 
Thus attempts at frequency analysis were abandoned, 
b. Randomly Rough Plane Surface 

As reported by Bailie [Ref. 34] and Hollis [Ref. 

35] , the generation of a boundary wave as predicted by Equa- 
tion (24) is very sensitive to the height of the source rela- 
tive to the plane of the roughness elements; although less 
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so for the relative height of the receiver. A series of 
acoustic measurements were made starting with the source 
housing height at the upper limit of the mechanism. This 
corresponded to a reading on the micrometer dial of 16.84 
relative cm. The height of the housing was then decreased 
in incremental steps of 1 mm until the top of the hole in the 
source cap was completely below the gravel. The maximum 
amplitude at each height position is compared to find the 
height at which the boundary wave was apparently strongest. 

This corresponded to a micrometer reading of 15.84 cm. With 
the source set at this height, the receiver was adjusted in 
the same manner. It seemed, however, that as long as the 
receiver was below the tops of the gravel in the immediate 
vicinity but above the midpoint of the gravel; the received 
signal was insensitive to receiver adjustment. Therefore, 
with the source and receiver set as specified, data collection 
began and the heights were not changed regardless of the 
variation of the gravel at any particular location. 

Each run consisted of the average of 2000 pings 
at a particular position on the randomly rough plane surface. 
The received signal was amplified, then filtered by 2 cascaded 
bandpass filters set with a passband of from 500 Hz to 110 kHz. 
The upper limit of this passband was established at this value 
as a consequence of the sampling frequency being set at 320 
kHz. Thus, because of the Nyquist criterion it was felt that 
reliable data up to at least 100 kHz could be obtained without 
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aliasing . The filtered signal was then amplified again to 
boost the signal up to a level that was well within the opera- 
ting range of the A/D converter. A complete table of equip- 
ment settings is contained in Table III, Appendix A. Each 
run was transformed to the frequency domain using Fast Fourier 
Transform (FFT) techniques and stored on magnetic disk for 
later use. 

The use of the probe tube limited the amount of 
useful data contained in each signal. Since the probe tube 
was 7.8 cm in length, the first evidence of the standing wave 
that was generated within the tube occurred at approximately 
450 microseconds after the beginning of the receiver's response 
to the acoustic signal. This left only 145 usable samples at 
a sampling frequency of 320 kHz. Since the algorithm for the 
FFT used for processing required the use of the number of 
data points which identically equal to a power of 2, the number 
of data points used for the FFT was 128. This established the 
frequency resolution, Af, as 



if = sampling fregu e n ^ „ ^20 000 ^ 5 kHz (32) 

number of samples 128 

Once an arbitrary radial line was designated as 
the zero angle, data were taken at the 40 cm and 20 cm range 
positions for angles 

9 = 9° + n 30° n = 0,1,2,. ..,11 (33) 
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At the 30 cm and 10 cm range positions, data were taken for 
angles 



e = 15 ° + n 30 ° n = 0 , 1 , 2,. ..,11 ( 34 ) 



Once these rough plate data were obtained, refer- 
ence data were collected using a smooth aluminum plate with 
the receiver and the bottom, of the hole in the source cap 
flush with the surface of the plate. Reference data were 
taken at 10 cm, 20 cm, 30 cm, and 40 cm. These reference 
data constituted the spherically spreading volume wave needed 
to calculate the boundary wave ratio as expressed in Equation 
(24) . 

To calculate the boundary wave amplitude, both 
rough plate and smooth plate data were converted into rec- 
tangular components for each frequency (magnitude and phase) 
at each range. Thus, for a particular range and frequency 
on the rough plate 

«R = “r + “r <35) 



And for the same range and frequency on the smooth plate 

*s ' ' “s'" 



The boundary wave is computed at that frequency and range 

= <“r-“s> ■" 3(3r-3s) Sb = V 
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where 



1/2 

Ag = t( = R -“s> <®R 

and 

>|)g = tan 

Therefore, the ratio of the rough plate frequency amplitude 
to the smooth plate frequency amplitude at a particular range 
is 



(aR -ag) 



1/2 



^BWA 

^WA 



lip 






B 










(38) 



This is the value to be compared with the predictions as 
outlined in Equation (24). This ratio was calculated for 
each range, angle, and frequency and stored on magnetic 
disk for comparison with theory. 
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V. EXPERIMENTAL RESULTS AND ANALYSIS 



A. SEAMOUNT 

Figures 19 through 21 are plots of pressure versus time 
for various configurations of seamount/ocean modeling. The 
top graph of Figure 19 is the free field pressure issued by 
the source, with a duration of approximately 0.15 milliseconds. 

The field diffracted over the 2D (PE Contour) model is 
shown in the middle and is weaker than the free field pressure 
has about the same duration. The pressure diffracted by the 
3D seamount model is somewhat greater than for the contour 
wedge and its shorter duration may be due to destructive 
interference for the longer paths around the seamount. 

The 2D contour model is of interest because it is the 
physical representation of what the two dimensional parabolic 
equation assumes in its solution of the problem. Figure 23 
shows this approximation to the real world without ocean sur- 
face (top), with smooth ocean surface (middle), and with rough 
ocean surface (bottom) . At very short times, up to 50 micro- 
seconds, the field is pure diffraction travelling by the 
shortest path; from 50 microseconds to approximately 200 
microseconds analysis reveals that the scatter is by a single 
reflection from the ocean surface, a situation that does not 
exist in the real ocean experiment and which is therefore 
not studied. The analysis of the forward propagated field 
starts at 203 microseconds and goes on for 400 microseconds. 
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FREE FIELD PRESSURE 
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FIGURE 19. rrcGSure versus time for free field, 2D (PE contour), 
and seamount with no ocean 



PRESSURE DIFFRACTED OVER 2D (P.E. CONTQURJ MODEL 

^ 203 O. 3 7a(««n t«ci 
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FIGURE 20, Pressure versus time for 2D (PE contour) model, diffracted 



The ocean equivalent time range is from 0.37 sec to 1.09 
sec. Chapman and Ebbeson [Ref. 40] have stated that this is 
the time between the first and second pulse observed at sea. 
What is clear from Figure 20 is that a rough ocean causes 
a significant reduction in this pulse, compared to the smooth 
ocean situation. This proves; if the ocean surface plays a 
part in the shadowing by a seamount, that it must be an 
appropriately rough surface that is assumed in the model. 

In Figure 21 the three dimensional model of the seamount 
is considered. Again, it is the period of this signal from 
203 microseconds to 603 microseconds with the purely dif- 
fracted signal seen in the top graph, that is of interest. 

The smooth ocean again causes a very large reverberant signal 
compared to the diffracted signal. But this is unrealistic, 
as the bottom graph shows for the rough ocean surface. The 
signal for the rough ocean is significantly lower than for 
the purely diffracted signal. Since the sum of the squares 
of the pressure amplitudes for each time increment for which 
the signal is above background noise is proportional to the 
energy contained in the signal for that time, a comparison 
of these sxxms of squares for the purely diffracted signal in 
the top graph with the signal contained between 203 micro- 
seconds and 603 microseconds in the bottom graph gives a 
quantitative feel for the relative amounts of energy involved 
For the purely diffracted signal the sum of the squares of 
the pressure amplitudes was 0.58; whereas, the sura of the 
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FIGURE 21. Pressure versus time for scale model seamount, diffracted 
smooth ocean, typical rough ocean 



squares for the rough ocean pressure amplitudes was 0.34. 

Thus, the diffracted signal contains more energy than the 
signal which interacts with the rough ocean, for the 35 knot 
wind condition. 

In spite of the difficulties with the variable source, 
there was one portion of the frequency domain analysis which 
was apparently valid. This was because data concerning the 
use of either of the seamount models with a model ocean tended 
to have the same kinds of variability. Thus when a ratio of 
two frequency spectra involving an ocean is formed, the varia- 
bility tended to divide out. This is shown in Figure 22. 

Here a comparison is made between the 2D and 3D seamount 
models with different types of oceans over each. The dashed 
line reaches a fairly stable relative value of -3 dB . This 
was calculated by subtracting the frequency data obtained for 
the 3D model with the average rough ocean from the 2D model 
with a smooth ocean. Thus the decibels are relative. 

The solid line reaches a fairly stable level of +7 dB. 

This was calculated by subtracting the frequency data for the 
2D model with the average rough ocean from the 3D model with 
the average rough ocean. Here, again, the decibels are rela- 
tive. Thus in a pseudo equation approach 

dashed line 2D, Smooth Ocean - 3D, Rough Ocean = - 3 dB 

solid line 3D, Rough Ocean - 2D, Rough Ocean = + 7 dB 

2D, Smooth Ocean - 2D, Rough Ocean = + 4 dB 

(40) 
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RELATIVK SIGNAL CONTRIBUTIONS FROM OFF-AXIS ROUGHNESS 
EI.FMFNTS AND SCATTERING FROM ROUGH OCEAN 



The second equation implies that approximately +7 dB is 
being contributed by roughness elements on the surface of 
the 3D model which are off the aixs of transmission. With 
the 2D model there is no mechanism for acoustic energy that 
is not being transmitted directly along the acoustic axis 
to be scattered back into the transmission path; whereas 
with the 3D model there is, and it seems to be contributing 
a relatively significant amount of energy to the signal 
received in the shadow zone, near the seamount. 

The third equation suggests that in the far field the 
presence of the smooth ocean above the seamount increases the 
received signal by approximately +4 dB when compared to the 
same model with a rough ocean. 

B. RANDOMLY ROUGH PLANE SURFACE 

After the rough plate data and the reference data had 
been collected and recorded on magnetic disk, the boundary 
wave amplitude and phase and the ratio of the boundary wave 
to the reference wave amplitude and phase were calculated for 
each angle and each range for the frequencies 2.5 kHz, 5 kHz, 
7.5 kHz, 10 kHz, 12.5 kHz, and 15 kHz. Then all angles for 
each range and each frequency were averaged using N-1 weight- 
ing to give a mean and standard deviation for each frequency 
and each range. These data are contained in Table IV, Appendix 
B. 

With this processed data available there were four essen- 
tial comparisons that needed to be made with Tolstoy's theory 
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for hemispherical scattering elements in an effort to estab- 
lish if the theory could be extended to the case of random 
roughness. The first was to determine the meaning of kh ^ 1 
when this was reevaluated using the spatial correlation length, 
L, of the gravel instead of the center to center spacing, h, 
of the hemispherical bosses . The second was to investigate 
the power law governing the amplitude ratio as a function of 
frequency . 

The third was to investigate the power law governing the 
amplitude ratio as a function of range. The fourth was to 
attempt an empirical determination of the scattering parameter, 
E, of the randomly rough surface by comparing the data with 
the scattering parameter calculated for the hemispherical 
bosses . 

In the work done by Bailie [Ref. 37], Medwin et al . , 

[Ref. 38], and Hollis [Ref. 39] using the hemispherical bosses, 

-5 

the scattering parameter was e = ah = 8 . 88 x 10 m where h 
is the center to center spacing of the bosses. 

To extend this to the case of random roughness, the 
correlation length, L, was substituted for h. Therefore 

kL 1 (41) 

As can be seen in Figures 23 through 26, a sharp peak is 
reached at 10 kHz. If it is assumed that this frequency 
marks the limit of coherence 



71 



— (10 kHz) L 



1 



345 



(42) 



(2tt) (10 kHz) 



5 . 5 mm 



This is in good agreement with the measured correlation 
length of the gravel surface of 



Thus, the use of the spatial correlation length as defining 
the limit of coherent scatter in extending Tolstoy's boun- 
dary wave theory is appropriate. 

The next question addressed was the power law of the 
amplitude ratio as a function of frequency. The values of 
the amplitude ratio versus frequency on a log-log scale are 
contained in Figures 23 through 26. Using linear regression 
techniques, the slope of the best straight line fit from 2.5 
kHz to 10 kHz was calculated for each of the four ranges . 

The average slope for the amplitude ratio as a function of 
frequency was then calculated from these four slopes using 
N-1 weighting. That average slope was +1.50 which is in 
agreement with the predicted slope of 3/2. 

The third question addressed was the power law associated 
with the amplitude ratio as a function of range. Figures 27 
through 30 are the plots of amplitude ratio versus range on 
a log-log scale for each frequency from 2.5 kHz to 10 kHz. 



L 



4 . 8 mm 



(43) 
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FIGURE 23 
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FIGURE 24 
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figure 25 
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FIGURE 26 
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FIGURE 27 



77 



Bwa/Vwa 




OF RANGE. THEORY LINE USING SCATTERING 
PARAMETER FOR HEMISPHERICAL BOSSES h = 0.2 CM 



FIGURE 28 
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FIGURE 29 
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FIGURE 30 
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As can be seen from these figures, there is no clear agree- 
ment with theory in this case; although the average of the 
four slopes is +0.43. Since the agreement was so striking 
between theory and data for the case of amplitude ratio as 
a function of frequency, this cannot be discounted. More 
experimentation would be necessary before anything defini- 
tive could be said about the range dependence. 

Finally the question of an empirically derived scattering 
parameter was addressed. It was assumed that the average 
data as listed in Table III, Appendix B was the value of 
Tolstoy's theory for a particular range, frequency, and 
unknown scattering parameter 



[ 



BWA, 

VWA^ 



■R 



(2-nr) 



1/2 ^ 3/2 



(44) 



where 



r = a specified range 

k = wave number at a specified frequency 
£„ = unknown scattering parameter 

Then, if Equation (24) were evaluated at that same range and 
frequency then divided into the data value in question. 



Data Value 

Calculated Value 



£j^(2TTr) 



1/2 



e^(2irr) 



1/2 



k ^/2 




(45) 
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where 



e 



T 



8.88 ^10 ^ m 



as calculated for the case of hemispherical bosses with 0.2 
cm center to center spacing. The resultant number is then 
multiplied by the theoretical scattering parameter for the 
hemispherical bosses 



Calculated Value ^ 

This procedure was followed and a numerical value for the 
scattering parameter was found for each range and frequency. 
The resultant 16 numbers were averaged' using N-1 weighting. 
As reported in Ref. 16, the scattering parameter for the 
hemispherical bosses was 



where a was the rms height of the hemispherical bosses. 
Following this approach using the calculated statistics 
for the gravel 



Data Value 



0.300 a 



(47) 



R 



n a 



(48) 



where 



R 



average scattering parameter for 
gravel = 6.03 x 10~^ m 



n 



numerical coefficient which was 0.300 for 
the case of hemispherical bosses 
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then 



n = — = 0.27 (49) 

a 

This result implies that the scattering parameter as pre- 
dicted by Tolstoy for hemispherical bosses is the same for 
the randomly rough plane surface. 
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VI. CONCLUSIONS 



A. SEAMOUNT 

If the pressure versus time graphs in Figure 25 are com- 
pared, the first pulse observed in the 79 km shadowing sit- 
uation at sea [Ref, 40] was the direct diffracted signal and 
the second pulse was the surface scattered energy. The time 
separation observed in the laboratory between these two 
groups of energy correspond to that which Chapman [Ref. 41] 
and Ebbeson observed at sea. This concept is represented 
graphically in Figure 31. In this figure the time from the 
onset of the signal until 203 microseconds is from the graph 
of the energy diffracted over the 3D seamount model; while 
the time from 203 microseconds until 603 microseconds is 
from the graph of the 3D seamount model with a rough ocean. 
Second, if the entire graph of pressure versus time for the 
seamount with a smooth ocean is compared with that of the 
seamount with a rough ocean (the bottom two graphs of Figure 
24), the distinctly different nature of a reflection dominated 
situation and a diffraction dominated situation becomes m.ore 
apparent. When the ocean is smooth, multiple reflections are 
the principal source of energy. When the ocean is roughened 
by a 35 knot wind the diffracted pulse (the first pulse) is 
stronger than the multi-scattered, surf ace-to-seamount , 
energy (second and later pulses) . A CW experiment would be 
diffraction dominated with these seas. 
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FIGURE 31. Composite pressure versus time showing diffraction 

only first pulse and scattered second and later pulses 



B. RANDOMLY ROUGH PLANE SURFACE 

The work done with the randomly rough plane surface 
constitutes a first look at the problem of extending Tolstoy's 
boundary wave theory to the real world. The work done in 
this thesis suggests that Tolstoy's theory can be applied to 
the real problem of a randomly rough rigid surface. This is 
of fundamental importance because of the decreased loss by 
geometric spreading that the boundary wave suffers relative 
to smooth surface spreading. 

To this point the seamount and the randomly rough plane 
surface have been discussed as if they were separate problems. 
In reality they are not. The Dickens Seamount in the Gulf 
of Alaska has small scale roughness on its surface which has 
not been represented on the 3D model. Because of the great 
depths of Dickens Seamount this small scale roughness has 
never been observed. However, as stated by Chapman [Ref. 42] 
the spectra for the shot data collected at sea show a curious 
interference pattern at a relatively stable frequency differ- 
ence whenever ray diagrams indicate the acoustic wave is 
interacting with the seamount. As pointed out by Hollis 
[Ref. 43] for the case of a 14 degree wedge and hemispheri- 
cal bosses, there was a pronounced boundary wave diffracted 
over the crest of the wedge in his experiments. The phase 
speed for the boundary wave is slightly less than that of the 
volume wave [Refs. 44, 45]. It has been suggested by Medwin 
[Ref. 46] that the presence of the interference pattern in the 
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frequency data from DREP is evidence that the effects of 
the boundary wave may have already been observed in nature. 
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APPENDIX A 



TABLE II 

EQUIPMENT SETTINGS FOR SEAMOUNT EXPERIMENT 



A. Timing 

Word 

0 

1 

2 

3 

4 



Simulator 

Output 

1 

0 

0 

3 

0 



Time 

100 usee 
4 msec 
10 usee 
900 usee 
20 msec 



B. ARB 

1) Mode = trig 

2) block rate = 1 kHz 

3) fxinc = RAM 8 

X Y 

0 127 

1 -127 

0-255 0 

4 ) amp = 0.5V 

C. Polarization Voltage = +150 VDC 

D. Power Amplifier 



Q 


A 




V 


108 


2.5 




265 


PARS 








1) 


Gain = 


100 




2) 


Gain = 


20 





F. BP Filter 

1) Pass Band = 6 kHz - 40 kHz 

2) Pass Band = 20 Hz - 200 kHz 
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TABLE III 



EQUIPMENT SETTINGS FOR RANDOMLY ROUGH PLANE SURFACE EXPERIMENT 



A. ARB 

1) mode = trig 

2) block rate = 2.5 kHz 

3) func = tri 

4) starting addr = 126 

5 ) amp = 0.35 V 

B. Polarization Voltage = +150 VDC 

C. Power Amplifier = 220 V 

D. BP Filter 

1) Pass Band = 500 Hz to 110 kHz 

2) Pass Band = 500 Hz to 110 kHz 

E. Range Dependent Settings: 

R = 40 cm 

Timing Simulator 

word output time 



0 

1 

2 

3 

4 



1 


100 


ysec 


0 


1 


msec 


0 


400 


ysec 


3 


5 


msec 


0 


25 


msec 



R = 30 cm 



0 

1 

2 

3 

4 



1 

0 

0 

3 

0 



100 ysec 
1 msec 
120 usee 
5 msec 
25 msec 



PARS 



1 2 
50 100 



50 100 
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CjJNJMO LJtOMO 



R 



20 cm 



1 

0 

3 

0 



= 10 cm 



1 

0 

3 

0 



100 ysec 
800 ysec 
5 msec 
25 msec 



loo ysec 
510 ysec 
5 msec 
25 msec 



20 100 



20 100 
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APPENDIX B 



TABLE IV 

AVERAGE AMPLITUDE DATA BY RANGE AND FREQUENCY 



R = 40 cm 



R = 30 cm 



R = 20 cm 



R = 10 cm 



Freq (kHz) 


BWA (average) 


BWA/VWA 1 


2.5 


13.76 


0.41 


5 


11.04 


0.45 


7.5 


52.25 


1.84 


10 


35.86 


1.94 


12.5 


38.40 


1.23 


15 


16.88 


0.91 


2.5 


5.12 


0.15 


5 


13.81 


0.48 


7.5 


36.17 


1. 12 


10 


40.18 


1.69 


12.5 


43.92 


1.17 


15 


23.56 


0.92 


2 . 5 


5.01 


0.25 


5 


10.58 


0.59 


7.5 


27.74 


1.39 


10 


33.98 


1.56 


12.5 


32.06 


1.25 


15 


17.77 


1.16 


2.5 


6.36 


0.16 


5 


17.98 


0.51 


7 . 5 


36.20 


0.90 


10 


60.14 


1.36 


12 . 5 


59.46 


1.12 


15 


39.88 


1.24 
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